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Abstract 

Wild poliovirus type 3 (WPV3) has not been seen anywhere since the last case of 
WPV3-associated paralysis in Nigeria in November 2012. At the time of writing, the most 
recent case of wild poliovirus type 1 (WPVl) in Nigeria occurred in July 2014, and WPVl 
has not been seen in Africa since a case in Somalia in August 2014. No cases associated 
with circulating vaccine-derived type 2 poliovirus (cVDPV2) have been detected in Nigeria 
since November 2014. Has WPVl been eliminated from Africa? Has WPV3 been 
eradicated globally? Has Nigeria interrupted cVDPV2 transmission? These questions are 
difficult because polio surveillance is based on paralysis and paralysis only occurs in a 
small fraction of infections. 

This report provides estimates for the probabilities of poliovirus elimination in Nigeria 
given available data as of March 31, 2015. It is based on a model of disease transmission 
that is built from historical polio incidence rates and is designed to represent the 
uncertainties in transmission dynamics and poliovirus detection that are fundamental to 
interpreting long time periods without cases. 

The model estimates that, as of March 31, 2015, the probability of WPVl elimination 
in Nigeria is 84%, and that if WPVl has not been eliminated, a new case will be detected 
with 99% probability by the end of 2015. The probability of WPV3 elimination (and thus 
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global eradication) is > 99%. However, it is unlikely that the ongoing transmission of 
cVDPV2 has been interrupted; the probability of cVDPV2 elimination rises to 83% if no 
new cases are detected by April 2016. 

Added July 10, 2015: On June 26, 2015, a paralytic polio case caused by ongoing 
circulation of established cVDPV2 lineage was confirmed by the Global Polio Laboratory 
Network. The date of paralysis onset was May 16, 2015, 181 days after the previous case. 
The case provides new information about the prevalence of cVDPV2 in Nigeria that can 
be used to both assess the accuracy of the predictions for elimination and case detection in 
the first version of the paper and provide the basis for an updated forecast for the time to 
cVDPV2 elimination if no cases are detected. 

The model predictions for cVDPV2 pepared on March 31, 2015 proved to be accurate: 
elimination did not occur and the observed case occurred within the 95% prediction 
interval. The updated forecasts predict that this one case has little impact on the 
estimates for the probability of elimination by April 2016. This counter-intuitive result 
was found because the 6-month silent interval following the previous 13-day interval 
indicates that cVDPV2 prevalence has fallen roughly 16-fold in Nigeria since November 
2014 and that the recent mean effective reproductive number for ongoing cVDPV2 
transmission in Nigeria is below 1. Under the most likely scenario for the transmission 
dynamics in the near future, the median time to the date of onset of the next case if 
elimination does not occur is July 15, 2015. 

I intend to provide further updates to this manuscript if new cases occur. It is 
important to note that if another case occurs before the end of July 2015, it will likely 
indicate increasing prevalence and so forecasts for elimination by April 2016 would then be 
more pessimistic. 


Introduction 

Since the World Health Assembly announced the intention to eradicate polio worldwide in 
1988, all countries except Pakistan, Afghanistan, and Nigeria have interrupted indigenous 
transmission of all three wild poliovirus (WPV) serotypes at least once [T] . Nigeria may be 
the next country to achieve wild-type polio elimination. The most recent case of wild-type 
1 (WPVl) in Nigeria occurred in July 2014 [5]. The most recent case of wild-type 3 
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(WPV3) seen globally occurred in Nigeria in November 2012 [I]. Is Nigeria wild-polio-free? 
Is WPV3 eradicated globally? Furthermore, the last case of WPVl in Africa occurred in 
Somalia at the tail end of the Horn of Africa Outbreak in August 2014. Is Africa 
wild-polio-free? 

The primary mechanism for detecting poliovirus in a population is surveillance for 
paralytic poliomyelitis. The incidence of poliomyelitis is tracked by the acute flaccid 
paralysis (AFP) surveillance and global polio laboratory network [31|3]. The certification 
period to declare that a polio serotype has been eliminated from a previously endemic 
region with certification-standard surveillance is three or more years without a polio case 
associated with indigenous WPV circulation This criterion is based on experience with 
regional elimination [5] and is necessarily long because paralysis due to polio is uncommon. 
Estimated case-to-infection ratios in fully-susceptible individuals are 1:200 for wild type 1, 
1:1900 for wild type 2, and 1:1150 for wild type 3 [7]. Previous modeling work has 
supported that the three year rule of thumb is reasonable, but that more precise estimates 
of the elimination time depend on the case-to-infection ratio for each serotype and the 
local conditions leading up to elimination [HHin]- 

Polio elimination in the most challenging settings has been driven by supplemental 
immunization activities that require extensive resources and prolonged community 
engagement to be successful [UIIIIIII]. Operational capacity and community engagement 
can be difficult to maintain for years in the presence of political instability, humanitarian 
crises, and substantial health care needs unrelated to polio [HITS]. While it is reasonable to 
maintain the three year certification period as the gold standard, the operational, 
humanitarian, and financial requirements support the need for specific, data-driven 
estimates of the likely period of silent polio persistence. 

The need for specific silent duration estimates is especially important in the context of 
the plan to stop all use of trivalent oral polio vaccine (tOPV) in April 2016 and replace it 
with bivalent OPV (bOPV) containing only serotypes 1 and 3 to prevent the seeding of 
new type 2 circulating vaccine-derived polio (cVDPV2) outbreaks [THITS] . cVDPV 
outbreaks can occur when vaccine strains transmit for sufficient time to allow for genetic 
reversion of the markers of attenuation m- OPV2 is the most common source of cVDPV 
globally m, and Nigeria has had ongoing transmission of cVDPV2 since 2005 min]. In 
2014, cVDPV2 was the dominant cause of poliomyelitis in Nigeria [2], and this is likely 
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because Nigeria has been primarily using bOPV to drive the successes in WPV 
control [TSlfT^ . However, due to expanded use of tOPV and combined tOPV/inactivated 
polio vaccine campaigns from the second half of 2014 muo], cVDPV2 cases in Nigeria 
abruptly stopped appearing at the end of November 2014 after producing 30 cases earlier 
that year [2]. Are the established cVDPV2 lineages gone? Will they be gone before the 
planned April 2016 tOPV cessation deadline? When is it reasonable for Nigeria to switch 
from a strategy focused on WPV elimination to one focused on cVDPV2? 

To estimate the probabilities of elimination from Nigeria given the data available at the 
time of preparation, March 31, 2015, 1 built a model appropriate to answer the question; 
how long do we have to wait before a chain of transmission either terminates or produces 
another case? The model is based on four premises: (1) that elimination only depends on 
the last few hundred to few thousand infections in a much larger population, (2) that we 
only need to see one more case to know elimination did not occur, (3) that we do not care 
where infected people are located because AFP surveillance will find any cases in 
proportion to the surveillance sensitivity, and (4) that the few remaining pockets of 
infection have similar transmission dynamics to the many sources of polio cases in the past. 
The model incorporates both the parametric and stochastic uncertainties that limit our 
knowledge of polio transmission. It provides estimates of the probability of elimination 
given no new cases and the time to the next case if elimination does not occur, and the 
results are appropriate for supporting the rational analyses m of planned polio 
vaccination policies. 

Added July 10, 2015; On June 26, 2015, a paralytic polio case caused by ongoing 
circulation of established cVDPV2 lineage was confirmed by the Global Polio Laboratory 
Network. The date of paralysis onset was May 16, 2015, 181 days after the previous case. 
With this new information, I assess the accuracy of the predictions for elimination and 
case detection in the first version of the paper and provide an updated forecast for the 
time to cVDPV2 elimination if no cases are detected. New sections are marked by the 
“Added July 10, 2015” tag. 
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Results 


The left column of Fig. [T] depicts elimination: it shows the probability by date that each 
serotype been eliminated if no new cases are detected. The right column depicts case 
detection: it shows the probability of the time to the next case if elimination does not 
occur. The dashed line indicates March 31, 2015. The solid blue curves provide estimates 
based on models with mean effective reproductive number equal to 1 (mean growth rate of 
the number of active infections equal to 0). For a fixed surveillance sensitivity, Red = 1 is 
the worst-case scenario because it maximizes the possible duration of silent transmission 
without elimination. Deviations from this assumption either make elimination easier 
{Reff < 1) or make the silent period between cases shorter {Reff > !)■ For WPVl and 
WPV3, the green curves provide less conservative estimates based on the assumption that 
vaccination quality maintains mean i?eff < 1 (growth rate < 0) at the level estimated at 
the time of the last case. For cVDPV2, the green curve assumes that the mean growth 
rate is held at one standard deviation below zero. Comments below summarize results 
from the green curves unless otherwise stated. For WPVl and WPV3, the dashed blue 
curves provide conservative estimates under the assumptions that i?eff = 1 and surveillance 
is only likely to detect 50% of polio cases. The poor surveillance curves are not shown for 
cVDPV2 because elimination is very unlikely, occurring in less than 1 in 100 000 
simulations. 

WPVl. The most recent WPVl case in Nigeria occurred on July 24, 2014. As of March 
31, 2015, there is an 84% chance that WPVl has been eliminated from Nigeria. The chance 
of elimination will be 93% by the end of 2015 if no new cases are seen even under the 
worst-case scenario for surveillance and the force of infection (Fig. [T]4, dashed blue curve). 
If WPVl does not eliminate before making another case, there is a 66% chance we will see 
a new case by the end of June 2015 and a 95% chance by the end of 2015 (Fig. [TJ3). 

The most recent case outside of Nigeria occurred in Somalia on August 11, 2014 at the 
tail end of the Horn of Africa outbreak that begin in 2013, 69 days after the case preceding 
it [22]. The corresponding predictions from the model for Somalia are thus essentially 
identical to those for Nigeria. The combined probability that both countries have 
eliminated (and thus all of Africa) is the square of the Nigeria values, and so the model 
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Figure 1. Model results for elimation and time to next case. Probability of 
elimination assuming no new cases are observed: WPVl (A), WPV3 (C), cVDPV2 (E). 
Probability of observing a new case if elimination doesn’t occur, WPVl (B), WPV3 (D), 
cVDPV2 (F). Fig. prepared March 31, 2015 (black dashed line). In each panel, the 
horizontal axis origin is the time of the most recent case, the solid blue curve shows the 
scenario with mean iieff = 1 and perfect surveillance, the green curve shows a less 
conservative scenario with mean i?eff < 1 and perfect surveillance, and the dashed blue 
curve shows a conservative scenario with mean i?eff = 1 and 50% surveillance sensitivity. 
For cVDPV2 (C,F), the orange curve depicts an optimistic scenario in which the mean 
i?eff is held at the lowest value ever observed and the standard deviation is reduced to 
one-fourth of its observed value. 


estimates a 70% chance that Africa is wild-polio-free as of March 31, 2015. 

WPV3. The most recent WPV3 case globally occurred on November 11, 2012 in Nigeria. 
As of March 31, 2015, there is a > 99% chance that WPV3 has been eliminated from 
Nigeria. The most recent case detected anywhere outside of Nigeria was in Pakistan seven 
months earlier (April 2012). The detection of the other serotypes since 2012 in both (and 
neighboring) countries indicates that surveillance is able to detect WPV3 if it were present, 
and so elimination from Nigeria likely represents global eradication. In the unlikely event 
WPV3 has not been eliminated, there is a > 99% chance we will see a new case by the end 
of 2015. 
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cVDPV2. The most recent cVDPV2 case in Nigeria occurred on November 16, 2014 
from a genetic lineage that has been in circulation since 2005. As of March 31, 2015, it is 
unlikely that the known cVDPV2 lineages nin] have been eliminated from Nigeria. 

If the known cVDPV2 lineages are still in circulation, there is an 84% chance the next 
case will appear by the end of June 2015, and a 97% chance a case will appear by the end 
of 2015. This model prediction for silent persistence is consistent with the longest observed 
interval of approximately 7 months seen from November 2012 to June 2013 [niiiHiiis]. 

The prediction for the probability of cVDPV2 elimination is sensitive to assumptions 
about future type 2 immunity. The rapid disappearance of cases in late 2014 (from one 
case every 10 days on average in 2014 to no cases in over 4 months) implies that type 2 
immunity is higher in the relevant population than it has been at any time previously, and 
so models based on historical transmission rates are likely too conservative. As population 
immunity affects transmission dynamics through herd effects that reduce the force of 
infection [53], consider three scenarios for the force of infection in 2015 and beyond: 

• Historically typical tOPV coverage (blue): For the historically typical coverage 
scenario, the model assumes that the future force of infection is described by the 
historical estimates—the mean growth rate is zero (i?eff = 1) and the standard 
deviation is described in the Methods. Under this scenario, elimination is unlikely 
with a 1% chance prior to June 2015. The estimated probability of elimination 
without seeing a case rises to 66% by April 2016, but does not reach 95% until June 
2017. 

• Good tOPV coverage (green): For a good coverage scenario, the model assumes that 
tOPV coverage improvements will maintain a reduction in the mean force of 
infection to one standard deviation below the historical norm and so not at 
unprecedented levels, A — /r = —3.44yr“^ (such that mean Reff ~ 0.85). Under this 
scenario, there is only a 3% chance of elimination by the end of June 2015 and an 
83% chance that cVDPV2 will be eliminated by April 2016 if no new cases are seen. 
The probability of elimination rises to 98% by the end of 2017. 

• Excellent tOPV coverage (orange): For the excellent coverage scenario, the model 
assumes that tOPV coverage improvements will maintain the mean growth rate at 
the lowest level transiently observed in the last ten years, X — n = —9.21yr“^ (such 




that mean Reff ~ 0.60), and will also reduce the variability of the growth rate to 
one-fourth of the historical standard deviation. Under this scenario, there is a 13% 
chance of elimination by the end of June 2015 and a 99% chance cVDPV2 will be 
eliminated by April 2016 if no new cases are seen. 

This analysis is focused solely on estimating the elimination probability for established 
cVDPV2 lineages and does not account for the emergence of new cVDPV2 lineages. In the 
event of a new emergence, these estimates would remain valid for the current genetic 
lineages, but additional modeling would be required to produce estimates for any new 
lineages and thus cVDPV2 overall. 

Added July 10, 2015: The new cVDPV2 case on May 16, 2015 provides information 
to validate the previous model and to update the cVDPV2 forecast. The original forecast 
was accurate in that it both correctly predicted that elimination was very unlikely prior to 
May 2015 (Fig. [1]^^), and the observed case occurred at the 54*^-percentile of the 
prediction interval for the time to the next case given the data available as of March 31, 
2015 (Fig. [2^). 

For each of the three scenarios considered for the original forecast, I estimated the 
probablility that each scenario produced the observed outcome. For the historically typical 
tOPV coverage scenario (blue), the observed outcome of a case after 181 days occured in 8 
of 100 000 simulations; for the good tOPV coverage scenario (green), 23 of 100 000 
simulations; and for the optimistic scenario (orange), 4 of 100 000 simulations. The 
absolute probabilites are low because any specific point outcome is unlikely, but the 
differences in the probabilities provide evidence to weight how consistent each scenario is 
with the observed outcome. That the green scenario is most likely to produce the observed 
outcome supports the assumption of the original paper that a scenario in which the mean 
Ref f is below 1 but not at unprecedentedly low levels is reasonable to describe the present 
situation. 

As described in the Methods, 1 used the set of previous simulations consistent with the 
observed outcome to estimate the prevelance of cVDPV2 infection in Nigeria on May 16. 
The median estimated infection prevalence on May 16, 2015 is 103 infections, which is 
roughly 16-fold lower than the median estimated prevelance in November 2014. Thus, the 
new case, while unfortunate, provides evidence of significant reductions in cVDPV2 
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Figure 2. Comparison of observed May 16, 2015 cVDPV2 case with the 
original forecast, and updated forecasts. The dashed black line marks the date of 
onset of the most recent cVDPV2 case on May 16, 2015. (A) Original forecast from March 
31, 2015 for the probability of cVDPV2 elimination assuming no new cases are observed 
and (B) the probability of observing a new case if elimination doesn’t occur (reproduced 
from Fig. [T|). (C) Updated forecast prepared on June 30, 2015 for elimination if no new 
cases are observed and (D) probability of the time to the next case if elimination does not 
occur. (E) Overlay of new forecast (solid lines) and original forecast (dashed lines) for 
easier comparison. As in Fig. [U the solid blue curve shows the conservative scenario with 
mean i?eff = 1 and perfect surveillance, the green curve shows the scenario most consistent 
with the recent 6-month silent period with mean i?eff < 1 and perfect surveillance, and the 
orange curve depicts an optimistic scenario in which the mean i?eff is held at the lowest 
value ever observed and the standard deviation is reduced to one-fourth of its historical 
value. Starting from May 16, 31% of all simulations in the blue scenario eliminate without 
producing another case, 61% eliminate in the green scenario without a case, and 73% 
eliminate in the orange without a case. 


prevalence. 

From the estimated initial condition distribution, the model produces updated forecasts 
from May 16, 2015 for the probability of elimination given no new cases (Fig. [Sp) and the 
time to the next case if elimination does not occur (Fig. [2jl)). Under the most likely of the 
three simulated scenarios for Reff in which the mean effective reproductive number is one 
standard deviation below the historical norm (green), the probability of elimination if no 
new cases occurs is 64% by the end of 2016, 85% by April 2016, and 99% by the end of 
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2016. These results are essentially unchanged from the estimate prepared on March 31 
(Fig. Upl). The model actually increases its estimate of the probability of elimination in the 
next few months because the new case reduces the uncertainty in the space of possible 
trajectories that allow for persistance through May without cases or elimination. If 
elimination does not occur, the median predicted date of onset of the next case is July 15, 
and there is a 92% chance of another case by the end of the year (Fig. HP). 

Discussion 

As of March 31, 2014, the model predicts that WPV3 has been eradicated globally and 
there is an 70% chance that Africa is completely wild-polio free for the first time in history. 
Furthermore, in the now unlikely situation that WPVl has not yet been eliminated from 
Nigeria and that it persists with the worst-case mean reproductive number of exactly 1 
and only 50% surveillance sensitivity, we can expect to see it or for it to eliminate with 
93% certainty by the end of 2015. Thus, this author believes it will be reasonable to 
conclude that Africa is wild-polio free if no new WPV cases are seen by the end of 2015. 
This estimate for time to elimination of at most 14 months is significantly shorter than 
estimated in previous work [HHTQ] because the model takes into account the low infection 
prevalence at the time of the most recent case that is implied by the long interval between 
the most recently observed cases. 

For cVDPV2, despite the dramatic reduction in incidence implied by the lack of new 
cases since November 2014, we must be cautious about concluding too soon that the 
established lineages are gone. The low case to infection ratio of type 2 makes it the most 
difficult strain to detect at low prevalence, and the high rate of incidence in 2014 implies 
that the reservoir of existing silent infections in November 2014 was approximately 22 
times larger for cVDPV2 than for WPVl at the time of its last case. The larger reservoir 
delays elimination even with sufficient herd immunity. Furthermore, continued tOPV use 
may lead to new cVDPV emergences in populations with insufficient herd immunity. The 
disappearance of cases indicates that substantial improvement in type 2 immunity has 
significantly reduced the national incidence rate, but the historical stability of the effective 
reproductive number near 1 suggests that the dynamics in the populations that continue 
to transmit are more stable than the national trend overall. Estimates here apply only to 




the known endemic lineages and new emergences may change expectations for the time to 
complete cVDPV2 elimination even if no new cases are seen from the established lineages. 

The above points suggest to this author that it would be reasonable for the polio 
program in Nigeria to place more emphasis on cVDPV2 prevention by the end of 2015. 
However, until the world is wild-polio-free, the risk of the re-importation of WPVl must 
be balanced with the risk of ongoing cVDPV2 transmission and new cVDPV emergence. 
Furthermore, certification of polio-free status is a stringent goal with additional 
requirements for surveillance and specimen containment [5], and so these estimates need to 
be considered in that context. The model was designed to help policymakers manage that 
balance. Involved parties may use the probabilities over time to weight different scenarios 
and inform their planning m, and the model can be easily extended to consider 
alternative scenarios about the future dynamics of the force of infection. 

If a new case is detected, then these model predictions are no longer valid. However, 
since the silent periods as of March 31, 2015 are much longer than typical, a new case 
would likely represent one of two scenarios. It would either be the last orphan case very 
near elimination, or it will be the first case of a new outbreak for which more cases will 
shortly follow. In the event another case occurs but is not followed by an outbreak, this 
analysis can be repeated to provide updated estimates. In the event of an outbreak, 
surveillance sensitivity is not the primary policy concern. 

The model assumes that the balance of birth and vaccination maintains stationary 
dynamics for the force of infection with constant mean infection growth rate and random 
variability. With additional data about birth demographics and vaccination activities, this 
model can be extended to include non-random structure in the statistical model of the 
force of infection. One could inform such models by correlating the apparently random 
changes in the force of infection with the supplemental vaccination calendar [2l [Tn[T^f20] 
to estimate how future campaigns may affect the force of infection and subsequent 
probability of elimination. 

The model is based solely on AFP surveillance and ignores any role for alternative 
active surveillance strategies such as environmental surveillance [25j . Certification-quality 
AFP surveillance provides a representative sample from all populations regardless of 
location. It is not obvious how to include environmental surveillance (ES) in the model 
because ES is only sensitive at fixed locations, but the model makes no assumptions about 




spatial structure. Furthermore, at a given ES site, there is not yet any practical way to 
quantify the meaning of a negative sample. And while positive environmental samples 
provide unambiguous evidence of local transmission at the time of the sample, we do not 
yet know how to quantify the effect of a positive sample on quantitative estimates of 
prevalence and thus future elimination probabilty. Because ES can dramatically improve 
the sensitivity in well-characterized populations [ISJUn] > the role of ES in certifying 
elimination warrants further empirical and theoretical study. 

Added July 10, 2015: After this paper was submitted, the GPEI reported an 
environmental sample from March (Kaduna State, Nigeria) that was positive for cVDPV2 
from an established lineage and a new case from the same lineage with a date of onset 
of Mayl6, 2015 was reported on June 30 in the Federal Capital Territory, confirming the 
predictions of this paper that cVDPV2 elimination was unlikely at that time and that a 
new case was likely within a few months of March 31, 2015. 

The model estimates that this one case does not significantly change the forecast for 
cVDPV2 elimination by April 2016. This perhaps counter-intuitive result occurs because 
the new case is consistent with a substantial drop in cVDVP2 prevalence over the last 6 
months. The median predicted time to the next case is 2 months under the most likely 
scenario for the future effective reproductive number. If another case with a date of onset 
prior to the end of July 2015 is observed, then it will be likely that cVDPV2 prevalence is 
rising, and this would lead to significantly more pessimisitc forecasts for cVDPV2 
elimination by April 2016. 1 will continue to provide updated forecasts if new cases occur. 


Methods 

The persistent intervals prior to elimination and the silent periods between cases given no 
elimination are determined by the force of infection, the case-to-infection ratio for each 
serotype, the mean infectious duration, and the number of extant infections at the time of 
the most recent case. In contrast to complex models with many elements that explicitly 
model assumptions about the influences of demographics, contact patterns, heterogeneous 
immune states, and vaccination history on transmission [SHin], this model reduces all the 
complexity of transmission into a statistical model for the time-varying force of infection. 
This design facilitates efficient marginalization over the substantial uncertainty in 


12/M 




transmission dynamics in a manner consistent with the limited available data. The range 
of the simulated outcomes from the model honestly represents the contributions of both 
parametric and stochastic uncertainties to the estimates of the probabilities of elimination. 

Dynamical model. The deterministic approximation to the model is given by a single 
equation: 

iV, = (A(t)-/r)7V, 

where is the number of infected people at time t, X(t) is the time-varying force of 
infection, and /i is the inverse of the mean infectious duration. The growth rate is 
(A(t) — fj,) and the effective reproductive number is i?eff = ^- In. th® stochastic version of 
the model, new infections occur at the inhomogeneous Poisson rate X(t)Ni{t) and existing 
infections are cleared at Poisson rate /i. 

This simple linear model is appropriate to model transmission near elimination because 
we can make a few assumptions that are not generally true. (1) There are many fewer 
infected than susceptible people and so the force of infection is independent of the 
susceptible fraction. (2) To know if elimination does not occur, we only need to predict the 
next paralytic case and do not need to correctly model outbreak size or duration. (3) Since 
we do not need to model the non-linear dynamics during an outbreak, but we want to 
explore the wide range possible transmission histories that could occur between cases or 
prior to elimination, a flexible statistical model of the force of infection is more 
appropriate than a specific dynamical model with mechanistic assumptions. 

Statistical model for the time-varying force of infection. A statistical model of 
the force of infection was constructed from published polio incidence data (ISl Table!) . The 
paralytic polio cases were re-binned to count the number of cases, C(t), in At = 3 month 
intervals and the growth rate between intervals was estimated as 

X{t + At) — M ^ log ■ When there are no cases in an interval, the growth rate 

was assumed to be the same as in the nearest bin for which there was data (Fig. EJ. 

The distribution of growth rates is approximately normal. There are no significant 
differences by serotype (Fig. |4]4) or significant deviations from normality (KS-test, 
p > 0.05). The mean growth rate is —0.38yr“^ and insignificantly different from 0; the 
standard deviation is 3.44 yr“^; there is no significant dependence of the growth rate with 
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Figure 3. Historical polio incidence in Nigeria. Polio case counts in 3 month bins: 
WPVl (A), WPV3 (C), cVDPV2 (E). Inferred growth rate (and reproductive number for 
/i = 26yr“^): WPVl (B), WPV3 (D), cVDPV2 (F). Note that while the overall case 
counts vary over time by two orders of magnitude, the reproductive number mean and 
range are stable. This motivates the assumption that transmission dynamics are stable 
among the populations that support ongoing transmission even as the number of such 
populations is lower now than a decade ago. 


time. It is more intuitive to consider the effective reproductive number. For a mean 
infectious duration of 2 weeks, fjL = 26yr“^, the estimated mean effective reproductive 
number is 0.99 as expected for an endemic disease, and the standard deviation is 0.14. To 
test for temporal correlations indicative of seasonality, I examined the Fourier power 
spectra of each growth rate time series. The spectrum for WPVl indicates the existence of 
seasonality with peak power at a 1 year period, but there are no clear peaks for WPV3 or 
cVDPV2 (data not shown). Despite decreases in polio prevalence and increases in 
vaccination quality since 2004 [28], the observed growth rate mean and range appear to be 
stable. 

The estimator is insensitive to changes in the surveillance quality over time as long as 
the fraction change in surveillance sensitivity, (5s, in the estimation time interval is small. 
The growth rate estimator with a term for the unknown changing surveillance sensitivity 
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Figure 4. Description of the model. (A) Quantile-quantile plot of the growth rate 
distributions for each serotype (Fig. [3]B,D,F) vs a normal distribution with mean 0 and 
standard deviation 3.59 yr“^. There are no significant differences between serotypes or 
deviations from normality (KS-test, p > 0.05). (B) Dynamical model cartoon: the model 
simulates the fade out and persistence of chains of transmission among a small number of 
people in a larger population. (C) Example model trajectories. Simulations were 
terminated when either enough infections were created to cause a new case or when 
elimination occurred. (D-G) Example simulated growth rate trajectories for the model 
(also indicated is i?eff for p = 26yr“^). Compare with Fig. [3]B,D,F. 


is: log ^ ^ ^ ^ S' typical size of the first term is 

the standard deviation of the growth rate, and so as long as Ss 0.86, the change in 
surveillance sensitivity can be ignored. The national non-polio AFP detection rate has 
roughly doubled since 2005 mm- Thus, while the slow increase in surveillance sensitivity 
is significant over the ten year period, its impact on the statistical model for the force of 
infection is negligible {6s « 0.03 per quarter). 

To model possible trajectories for the time-varying force of infection, I constructed time 
series for the growth rate X{t) — p from independent samples for each 3 month interval 
from the growth rate distribution. The assumption of independence between intervals is 
justified because weak seasonality has little influence on the interval between cases (Fig. 4 
of ref. i)- Example trajectories are shown in Fig. |4l The mean infectious duration was 
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drawn from a uniform prior for fi on the interval fj, = [(30days)“^, (7days)“^]. The range 
of the prior is determined from the range of observed mean poliovirus shedding durations 
in various historical and immunological settings [30) . 

This statistical model represents changes in the balance of vaccination coverage, 
demographic dynamics, and relevant heterogeneities as assumptions about the mean, 
variance, and temporal correlation of the force of infection. It is based on the assumption 
that there are no long-term trends in the force of infection. This assumption is 
conservative in the sense that trends in either direction shorten duration of silent 
transmission by either making outbreaks or elimination more likely. Finally, note that this 
model is designed to simulate disease dynamics near elimination. The model would not 
work well if the saturating nonlinearities that determine outbreak size and duration were 
important because the model is inherently linear. 

Model of the interval between cases. To model the interval between cases, 
simulations were initialized at the time of the most recent observed case. I assumed that 
progression from infection to paralysis can be modeled as a Bernoulli trial with probability 
p = (case-to-infection). Accordingly, the probability that infection number Ni is the next 
paralytic case is geometrically distributed, P{Ni]p) = p(l — p)^'~^ [21]. Estimated 
case-to-infection ratios for individuals fully susceptible to paralysis are 1:200 for WPVl, 
1:1150 for WPV3, and 1:1900 for cVDPV2 |7] under the assumption that cVDPV2 has an 
indistinguishable phenotype from wild type 2 |15j . The impact of imperfect AFP 
surveillance with only 50% sensitivity (such that half of all polio cases can be missed on 
average) was modeled by reducing the case-to-infection ratios to 1:400 for WPVl, 1:2300 
for WPV3, and 1:3800 for cVDPV2 . 

Since the case-to-infection ratios above describe paralysis-susceptible infections only, 
and the statistical model for the force of infection is based on paralytic cases only, the 
model of new infections explicitly represents only infections in the cohort of people who 
are susceptible to paralysis. It does not represent infections in individuals that are 
protected from paralysis by prior immunity but are not fully protected from infection |7]. 
The model is agnostic to the details of how transmission between paralysis-susceptible 
individuals is mediated, and the role of individuals who can be infected but cannot be 
observed is subsumed into the empirical model for the force of infection between 
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paralysis-susceptible individuals. 


Initial condition distribution. The interval between the two most recent observed 
cases provides a target to calibrate the initial prevalence. The range of initial infections 
was calibrated to match the median interval between cases in the model to the observed 
interval between the two most recent cases. This is based on the assumption that prior to 
observing a silent period, the next case interval is expected to be similar to the previous 
one. 

For WPVl, the most recent case occurred on July 24, 2014 and the observed interval to 
the case prior to it was Ate = 58 days earlier [32] ■ (For comparison, since the start of 2012, 
the typical interval between WPVl cases was 5.5 [1, 27] days (mean [95% interval]).) For 
WPV3, the most recent case occurred on November 11, 2012 and the interval to the 
previous was Ate = 120 days [32] (vs. 15 [1,42] days). For cVDPV2, the most recent case 
occurred on November 16, 2014 and Ate = 13 days [32] (vs. 22 [1,141] days). 

For each fi, gradient descent was used to reduce the absolute difference between the 
median time to the next case in the model and the target interval to at most one day. The 
resulting fits of the initial conditions to the median intervals are closely approximated by a 
simple curve that depends on /i, the case-to-infection ratio for each serotype, the target 
interval between cases, and a factor of order 1; 


Nm - f 


(case-to-infection) 
^Ate 


-1 


with /wpvi = 0.55, /wPV3 = 0.60, and /cvdpv2 = 0.90. 

The initial condition curves indicate that at the time of the most recent case, the 
prevalence of cVDPV2 was approximately 22 times greater than the prevalence of WPVl 
and 73 times greater than WPV3 at the times of their most recent cases. 

Added July 10, 2015: Bayesian data assimilation for updated cVDPV2 
initial conditions. For the updated cVDPV2 forecast, a Bayesian data assimilation 
procedure was used to identify appropriate initial conditions. From the set of simulations 
from the March forecast that are consistent with the data as described in the Results, I 





built a mixture distribution for the initial number of infections on May 16, 2015: 


piNmW)= ^ P(scenario) P(lVi(0)[scenario, fj) 


scenario 


with 


P(blue) = 


35 ’ 


23 


35 ’ 


P(scenario) = p(green) = 

P( orange) = 

and normally-distributed initial conditions given the scenario 


P{NM 


scenario 


io,^) = < 


(IV,(0)1 blue,/r) ^jv(6j§^,200), 
(A^,(0)[green, - A/'(o.55^, 65^ , 
P(iV,(0)[orange, Ai) ^ AT(^0.09^, 25) . 


Simulation procedure. The stochastic version of the model was coded in C# and 
trajectories were simulated using the Gillespie’s direct method for time-varying 
rates [33l|34]. For each run of the model, the parameters A(t), and the number of 
infections required to produce the next case were sampled from their distributions, and the 
simulation was stopped when elimination occurred or the number of infections required to 
produce a case was reached. For all scenarios, at least 20 000 simulations and up to 
200 000 simulations were run to produce Fig. [T] 


Probabilities of elimination and the time to the next case. The probability of 
elimination without seeing another case was estimated as: 

, , _ ne(t) 

^eliminationlCj U) ’ 

iLe,yi) T iipyi) 

where ne{t) is the number of simulations that eliminate at or before time t and np{t) is the 
number of simulations that persist past time t. The probability of the time to the next 
case given elimination does not occur and given the observed silent period through March 
31, 2015 is estimated as the number of simulations that produce a case after March 31, 
2015 but at or before time t over the number of simulations that produce a case at any 
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time after March 31, 2015. Added July 10, 2015: For the updated cVDPV2 forecast, 
The probability of the time to the next case given elimination does not occur is estimated 
as the number of simulations that produce a case after May 16, 2015 but at or before time 
t over the number of simulations that produce a case at any time after May 16, 2015. 


Supporting Information 

SI Table 

Polio case incidence data table for 2004—2014. The polio case incidence rates for for 
WPVl and WPV3 through September 2014 were scraped from the figures of 
refs. [niiiiiiMisH] with the aide of the DataThief software package (B. Tummers, 
DataThief III. 2006:http: //datathief . org/). The cVDPV2 incidence data was taken 
from the supplement of ref. through the end of 2011 and extended through September 
2014 with scraped data from refs. [T1I2H]. The remaining data after September 2014 was 
derived from the weekly tally sheets available at ref. [5^ . 
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